Universality of the Tangential Shape Exponent at the Facet Edge of a Crystal 
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Below the roughening temperature, the equilibrium crys- 
tal shape (ECS) is composed of both facets and a smoothly 
curved surface. As for the "normal" profile (perpendicular to 
the facet contour), the ECS has the exponent 3/2 which is 
characteristic of systems in the Gruber-Mullins-Pokrovsky- 
Talapov universality class. Quite recently, it was pointed 
out that the ECS have a "new" exponent 3 for "tangen- 
tial" profile. We first show that this behavior is universal 
because it is a direct consequence of the well-establised uni- 
versal form of the vicinal-surface free energy (p: surface gra- 
dient): f(p) = /(0) +jp+ Bp 3 + 0{p 4 ). Second, we give 
a universal relation between the amplitudes of the tangential 
and the normal profiles, in close connection with the univer- 
sal Gaussian curvature jump at the facet edge in systems with 
short-range inter-step interactions. Effects of the long-range 
interactions are briefly discussed. 

PACS numbers: 05.70.Jk, 64.60.-i, 64.60. Fr, 68.35.Md, 68.35.- 
p, 82.65.Dp 

A crystal surface with small average tilt relative to a 
crystal axis is called vicinal surface. Below the roughen- 
ing temperature Tr, the vicinal surface is well descibed 
by the terrace-step-kink (TSK) picture where the sur- 
face is made up of many "flat" terraces connected by 
steps. Since the step is a unidirectional linear object, 
the vicinal surface, as an assembly of steps, belongs to 
the Gruber-Mullins-Pokrovsky-Talapov (GMPT) univer- 
sality class. JjJ In this view, the vicinal-surface free en- 
ergy f(p) (per projected area) as a function of the step 
density p has the well-known form of the expansion, [jl]]|] 
f(p) = /(0) + dip + a 3 p 3 + 0(p 4 ). This form of surface 
free energy leads to, via the Wulff construction, |3) the 
universal exponent 3/2 of the equilibrium crystal shape 
(ECS, for short) near the facet edge, z ~ (Air) 3 / 2 , where 
we have chosen the z axis to be facet normal, and Ax to 
be the "normal" distance from a facet edge. Careful 
consideration of the crystal anisotropy gives us a univer- 
sal relation Q between the coefficients a% and 03. This 
relation leads to the universal Gaussian curvature jump 
at the facet edge PM|| whose physical origin is, there- 
fore, quite differentpl from the universal curvature jump 
at T R . § 

Very recently, Dahmen et al. JlCj] pointed out that, if 
we approach the facet edge "tangentially" (i.e., along the 



tangential direction at the facet edge), ECS behaves as 
z ~ (Ay) 3 where Ay is the distance from the facet edge 
along the tangential direction. They drew this conclu- 
sion from an exact calculation (essentially equivalent to 
Ref.8) on the body-centered-cubic solid-on-solid model 
]IT[ (BCSOS model, for short). It is an interesting prob- 
lem, then, to see whether this "new" exponent is uni- 
versal or not. Also, relations, if any, to the above-listed 
universal behaviors are worth to be explored. To clar- 
ify these points is the aim of the present paper. As will 
be seen below, the tangential exponent 3 is a direct con- 
sequence of the GMPT-form of expansion of the surface 
free energy, hence, is universal. Further, the universal re- 
lation between the expansion coefficients leads also to a 
universal relation between the amplitudes of the normal 
and the tangential ECS profiles. 

Consider a facet and its neighboring curved surface of 
the ECS with macroscopic size at a temperature T be- 
low Tr (of the facet, to be precise). We use the Carte- 
sian coordinates (x, y, z) with ^-direction chosen to be 
facet normal (downword, for convenience), which allows 
us to describe the ECS by an equation z — z(x,y). 
The surface gradient at a position {x, y, z(x, y)) is conve- 
niently expressed by the two-dimensional gradient vector 
p = (px,P v ) defined by 



P. 1 



dz 
dx 7 



Py 



dz 
dy' 



(1) 



By f(p) we denote the free energy per projected area |T^] 
of the surface with fixed mean gradient p. Introducing 
a "field" r\ conjugate to p, we can perform the Legendre 
transformation p — > r], f(p) — > f(rf), with 



f(v) = f(p) VP- 



(2) 



The Andreev free energy f(rj) directly gives us the ECS 
as 111 



(3) 



where r = (x,y) and A is a scale factor. In what fol- 
lows we consider "normalized" ECS by setting A = 1 in 
(^) . Then the field rj is just the two-dimensional position 
vector of the ECS. 
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Properties of the ECS near the facet edge are gov- 
erned by the small- |p| behavior of f(p) which has the 
GMPT-type form of expansion. Properly taking account 
the crystal anisotropy, we have the expansion of the form, 



Hp) = /(o) + i(0)\p\ + B(6)\ P \ 3 + o(| P | 4 ), 



(4) 



where we have introduced the angle variable 9 by 

p x = \p\ cos 8, p y = \p\sin9. (5) 

Physically, j(9) is the step tension and 9 the mean run- 
ning direction angle of steps on the vicinal surface with 
gradient p. Note that 9 is also the direction angle of the 
tangential line of the facet contour (Fig.l). In systems 
with only short-range inter-step interactions, the coeffi- 
cient B(9) is always positive and is universally given by 



B{9) 



n 2 (k B Tf 
6j(9) 



(6) 



where j(9) = 7(6*) + d 2 ~f(9)/d9 2 is the step stiffness. In 
association with the Legendre transformation, we have 
the following relation between the two-dimensional posi- 
tion vector r — (x,y) and the gradient p: 



df( P ) 
dp 



(7) 



Near the facet edge in the curved region, we obtain from 
(0) (neglecting 0(|p| 4 ) and higher-order terms) 



x 
V 



x c (9) 

vM - 



\p\ 2 [3B(9) cos 9- B'(9) sin 9] 
\p\ 2 [3B(9)sm9 + B'(9) cos0]. 



(8) 
(9) 



where the curve {(x c (9), y c (0))}e is the facet contour cor- 
responding to the zero-gradient limit of the curved re- 
gion. Explicitly, we have 



x c (9) = 7 (6») cos ( 
t/ c (0)=7(0) sin ^ 



■7' (9) sinf 
i(9) cost 



(10) 

(11) 



We should note that (|l^) and ( |Ti"| ) are precisely the 
equations determining the two-dimensional ECS (=facet 
shape) from 7(0) regarded as a one-dimensional interface 
tension. [jl3| 

Let us now discuss the ECS near the facet edge in 
detail. For a given ECS described by an equation z = 
z(x,y), we specify any point on the ECS surface by the 
two-dimensional position vector (x,y). For convenience, 
we take the xy plane (i.e., z = 0) to be the facet plane, 
which corresponds to putting /(0) = in (Q). Fix a 
point P on the facet contour, and choose the x and y 
axes so that the y-axis is parallel to the tangential line 
of the facet contour at P, and x-axis perpendicular to 
it (Fig. 2). With this choice of the coordinate system, 
we have 9 = at P and 7 '(0) = (see §) and (0)). 
Our task is to obtain the ECS profile close to the point 
P=(x c (0),y c (0)) = (x c (0),0). 



Along the x-axis 
from (||) and (^jj) 

\P\ = Px 



0, Ax 



x c (0)), we have 



(Ax 



,1/2 



(12) 



giving the "normal" profile with the well-known exponent 
9 X = 3/2: 



3v/3~B(0) 



(Ax 



.3/2 



(13) 



The profile ( |l3|) leads to the divergent behavior of the 
normal curvature n x « d 2 z/dx 2 ~ (Ax)~ 1//2 near the 
facet edge. In the light of the universal Gaussian cur- 
vature jump at the facet edge, Q the "tangential curva- 
ture" K y « d 2 z/dy 2 1 along the x axis, vanishes |6|,|l4| as 
(Ax) 1//2 [Gaussian curvature is a product of two pricipal 
curvatures, k x and Ky], We should note that a differ- 
ent exponent y for the tangential profile z ~ (Ay) 9y 
(Ay = y — y c (0)) with 8 y > 2 has already been implied 
by the vanishing of the tangential curvature at P. 

The actual value of 9 y can be derived as follows. Note 
that, along the tangential line (x — x c (0) = 0), 9 and \p\ 
are not independent but are constrained to satisfy, 
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7(0)6» 2 + 3B(0)|p| 2 =0, 



which can be derived by expandin. 
respect to 9 and \p\ (\p\ << 1 and 
the point P). Combining ( pT[ ) with 



(14) 

§) and @ with 
<< 1, very near 
dill), we obtain 



7(0) : 



\P\ = 



Ay 



V6 7 (0)B(0)' 



(15) 
(16) 



along the y-direction. Near P, —z = /(x c (0), Ay + y c (0)) 
is expanded to give 



^ = [7(0) + |7"(0)^ 2 ]|p| + i?(0)|p|= 
-A 2/ |p|0-x c (O)|p|(l-^ 2 ). 



Putting ((TH) and (|T|) into (|T|), we have 
2 = 1 (At,) 3 



(17) 



(18) 



giving 9 y = 3. We should remark here that the results 
( |l3| ) and ( |l8| ) apply to any point on the facet contour, 
because the point P has been chosen arbitrarily; hence, 
7(0) and 5(0) in (|l|) and (|l|) can be replaced by 7 (6») 
and B(9). 

We have another simple geometrical derivation of dl§| ) 
as follows. Very near the point P, we take a different 
point Q on the facet contour (Fig. 3). Along the x'- 
direction chosen normal to the facet contour at Q, ECS 
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profile has the form (|T^) with B(0) replaced by B(9) 
where 9 (\9\ « 1) is the tangent angle (relative to y- 
axis at P) of the facet contour at Q. Note that the facet 
contour near P is approximately a part of a circle with 
its radius being the curvature radius R which is propor- 
tional to the step stiffness 7(0). 15 For the normalized 
ECS we have R — 7(0), hence, by an elementary geom- 
etry we can relate Ax' (distance along the x'-direction) 
to Ay as 



Ax' 



(Ay) 2 {Ayf 



2R 



27(0) 



(19) 



Replacing this Ax' with Aa; in ([l3]) we reproduce fll8| ) 
[we can put B(9) ss B(0), for our purpose]. 

We thus have shown that the "new exponent" 9 y = 3 is 
a direct consequense of the well-established GMPT-type 
expansion of the vicinal surface free energy (^). Let us 
next discuss the "critical amplitudes" of the profiles. In 
the derivation of the universal Gaussian curvature jump 
at the facet edge, [p] the universal relation 



B{9) 



7r 2 (fc B r) s 
67(0) 



(20) 



which holds for any system (with only short-range inter- 
step interactions) is essential. This relation was originally 
derived in the coarse-grained TSK picture of the vici- 
nal surface. || and has been confirmed in several ways. 
^,|||Jl^,|lj For example, exact calculation for the BC- 
SOS model have verified ( pp[ ) for arbitrary 9. || Using 
(pp|), we obtain a universal relation between the ampli- 
tudes of the normal and tangential profiles as follows. By 
A x (9) and A y (9) we denote the amplitudes of the ECS 
profiles, namely, 



A x (9)(Ax) 3 / 2 (normal direction), 
A y (9)(Ay) 3 (tangential direction) 



(21) 



Restoring the scale factor A in (||), we have from (13) 
M) and (ET 



[Al(9)A y {9)]^ = \ 



37rfceT ' 



which means that at a fixed temperature, the quantity 
(A 2 , Ay) 1 / 3 is constant along the facet contour. The scale 
factor A can be determined, for example, from the mea- 
surable ratio k/g(9) 2 (oc X/(k B T) |j) where k is the 
curvature of the facet contour and a the scaled fluctua- 
tion width of a single step. |Q 

Before closing, we give a brief comment on the effect 
of the long-range inter-step interaction (~ 1/r 2 , with r 
being the ter-step distance) which has its origin mainly 
in the elastic deformation, and is important in discussing 
real crystal surfaces. |l^] It has been known that inclu- 
sion of g/r 2 interaction with positive coupling constant g 
does not modify the GMPT-type form of expansion, but 
merely renormalizes the coefficient B. Explicit form of 



the renomalized B has also been known. ||T^,^,|TJ Hence 
the exponent 9 y = 3 does not change with the g/r 2 inter- 
action. Since 7(0) is a quantity associated with a single 
isolated step and is not affected by the g/r 2 interaction, 
the relation between 7 and the renormalized B should be 
modified. If g is ^-independent (9: mean running direc- 
tion of steps) then the universal relation (^) still holds 
in a modified form. However, if g depends on 9, A x and 
A y will no longer be universally related to each other. In 
this case, ( p2[ ) may provide a way to determine the 9- 
dependence of g experimentally, by measuring the ratio 

[A 2 x (9)Aym 1/3 /[A 2 x (0)Aym 1/3 - 

To summarize, we have discussed the tangential pro- 
file of the equilibrium crystal shape near the facet edge 
below the roughening temperature. We have shown 
that the profile z ~ (Ay) 9y (Ay: distance along the 
tangential direction at the facet edge) with 9 y = 3 is 
a direct consequence of the Gruber-Mullins-Pokrovsky- 
Talapov type expansion of the surface free energy f(p) — 
/(O)+7(0)|p| +B(9)\p\ 3 + . . ., which implies that the ex- 
ponent 3 is universal. Further, we have presented a gen- 
eral relation between the amplitudes of the normal and 
the tangential profiles, which results from the known uni- 
versal relation between the coefficient B(9) and the step 
stiffness j(6) = 7(6*) + j"(d). 
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Figure Captions 

Fig.l: Facet and vicinal surface. Upper right is an 
"atomic scale" view of the surface near the facet edge 
in the curved region, which can be regarded as an assem- 
bly of steps forming a vicinal surface. 
Fig. 2: Choice of x- and y-axes at P on the facet contour. 
Fig. 3: Geometrical derivation of Eq.(18). For Q very 
near P, "normal distance" Ax' relates to "tangential dis- 
tance" Ay as Ax' = (Ay) 2 /2R, where R is the curvature 
radius at P. 
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Figure 3 Y. Akutsu, N. Akutsu and T. Yamamoto 



